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Three-dimensional foam flow resolved by fast X-ray tomographic microscopy 
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Thanks to ultra fast and high resolution X-ray tomography, we managed to capture the evolution 
of the local structure of the bubble network of a 3D foam flowing around a sphere. As for the 2D foam 
flow around a circular obstacle, we observed an axisymmetric velocity field with a recirculation zone, 
and indications of a negative wake downstream the obstacle. The bubble deformations, quantified 
by a shape tensor, are smaller than in 2D, due to a purely 3D feature: the azimuthal bubble shape 
variation. Moreover, we were able to detect plastic rearrangements, characterized by the neighbor¬ 
swapping of four bubbles. Their spatial structure suggest that rearrangements are triggered when 
films faces get smaller than a characteristic area. 


Foam rheology is an active research topic BHa , mo¬ 
tivated by applications in ore flotation, enhanced oil re¬ 
covery, food or cosmetics [5]. Because foams are opaque, 
imaging their flow in bulk at the bubble scale is challeng¬ 
ing. To bypass this difficulty, 2D flows of foams confined 
as a bubble monolayer, which structure is easy to visu¬ 
alize, have been studied. However, the friction induced 
by the confining plates may lead to specific effects [6], 
irrelevant for bulk rheology. In 3D, diffusive-wave spec¬ 
troscopy has been used to detect plastic rearrangements 
BE]. These events, called Tls, characterized in 2D by 
the neighbor swapping of four bubbles in contact, are of 
key importance for flow rheology, since their combination 
leads to the plastic flow of foams. Magnetic resonance 
imaging has also been used to measure the velocity field 
in 3D [9]. However, both these techniques resolve nei¬ 
ther the bubble shape, nor the network of liquid channels 
(Plateau borders, PBs) within a foam. In contrast. X-ray 
tomography renders well its local structure. However, the 
long acquisition time of a tomogram, over a minute until 
very recently, constituted its main limitation, allowing to 
study only slow coarsening processes [iQiin]. 

Here, we report the first quantitative study of a 3D 
foam flow around an obstacle. Such challenge was tackled 
thanks to a dedicated ultra fast and high resolution imag¬ 
ing set-up, recently developed at the TOMCAT beam line 
of the Swiss Light Source m- High resolution tomogram 
covering a volume of 4.8 x 4.8 x 5.6 mm^ with a voxel 
edge length of 5.3 /im could be acquired in around 0.5 
s, allowing to follow the evolving structure of the bub¬ 
bles and PB network. Our image analysis shows that 
the 3D foam flow around a sphere is qualitatively simi¬ 
lar to the 2D flow around a circular obstacle: we reveal 
an axisymmetric velocity field, with a recirculation zone 
around the sphere in the frame of the foam, and a neg¬ 
ative wake downstream the obstacle. Bubble deforma¬ 
tions are smaller (in the diametral plane along the mean 



FIG. 1. A plastic bead of 1.5 mm diameter glued to a 
capillary is placed in the middle of the cylindric chamber of 
22 mm diameter and 50 mm height. The acquired tomograms 
cover the central region with a volume of 4.8x4.8x5.6 mm^. 
A typical X-ray projection image is shown on the right. 


direction of the flow z) than for a 2D flow, thanks to 
the extra degree of freedom allowing an azimuthal defor¬ 
mation: bubbles appear oblate before, and prolate after, 
the obstacle. Finally, we were able to detect plastic rear¬ 
rangements, characterized by the neighbor-swapping of 
four bubbles and the exchange of two four-sided faces. 
Our observations suggest that those events are triggered 
when the bubble faces get smaller than a characteris¬ 
tic size around given by a cutoff length of the PB 
130 jam in the case of our foam. 

Experimental set-up - We prepared a foaming solution 
following the protocol of m- we mixed 6.6% of sodium 
lauryl ether sulfate (SLES) and 3.4% of cocamidopropyl 
betaine (CAPB) in mass in ultrapure water; we then dis¬ 
solved 0.4% in mass of myristic acid (MAc), by stirring 
and heating at 60°C for one hour, and we diluted 20 
times this solution. A few mL of solution was poured in 
the bottom of a cylindrical perspex chamber of diameter 
22 mm and height 50 mm, compatible with 180 degree 
tomography. Bubbling air through a needle immersed in 
this solution, a foam was created until it reached the top 
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of the chamber. The bottom of the cell contains a tube 
connected to the open air and closed by a tap. Control¬ 
ling the opening of the tap, we could obtain a slow steady 
flow of the liquid foam. Its mean velocity determined a 
posteriori by image analysis is equal to 'r’flow = 8 /im/s. 
While flowing, the foam is deformed due to the presence 
of an obstacle, a smooth plastic bead of diameter 1.5 mm, 
attached to a capillary to fix its position in the middle of 
the chamber (Fig. [^. 

The experiments were performed at the TOMCAT 
beamline of the Swiss Light Source. Filtered polychro¬ 
matic X-rays with mean energy of 30 keV were incident 
on a custom made flow cell (Fig. attached to the to¬ 
mography stage with three translational and a rotational 
degrees of freedom. The X-rays passing through the 
foam in the chamber were converted to visible light by a 
100 /im thick LuAGiCe and detected by a 12 bit CMOS 
camera. Typically 550 radiographic projections acquired 
with 1 ms exposure time at equidistant angular positions 
of the sample were reconstructed into a 3D volume of 
4.8 X 4.8 X 5.6 mm^ with isotropic voxel edge length of 
ps = 5.3 pm. Such a 3D snapshot of the flowing foam is 
acquired in tgcan = 0.55 s, ensuring that motion artifacts 
are absent since 4can < order to follow the 

structural changes of the foam during its flow around the 
obstacle, we recorded a tomogram every 35 seconds for 
approximately 20 minutes (resulting in around 36 tomo¬ 
grams). The tomograms quality is enhanced using not 
only the X-rays attenuation by the sample, but also the 
phase shift of the partially coherent X-ray beam as it in¬ 
teracts with the PBs and senses the electron density vari¬ 
ation in the sample m- This phase shift was retrieved 
using a single phase object approximation m- 

Image analysis - The tomograms are then segmented, 
separating the PBs and vertices from air. Fig. shows 
two successive time steps of the 3D snapshots of the re¬ 
constructed PBs network during the foam flow around 
the sphere. We measured the liquid fraction from the 
segmented images, by computing the relative surface oc¬ 
cupied by the PBs and vertices on individual horizontal 
slices. We measured an averaged liquid fraction of 4% 
over a tomogram, which did not evolve significantly dur¬ 
ing our experiments. 

Then, we reconstructed and identified individual bub¬ 
bles of the flowing foam, following the procedure we re¬ 
cently developed and validated on static foam samples, 
imaged at the same acquisition rate and spatial resolu¬ 
tion m- We did not observe any evolution of the size 
distribution of the polydisperse foam studied here, with 
an average volume V = 0.36 ±0.13 mm^, hence coarsen¬ 
ing remains negligible. Typically, 160 bubbles are tracked 
between two successive 3D snapshots, leading to statistics 
over 5600 bubbles. To ensure a high tracking efficiency, 
we took into account the mean motion of the foam. Bub¬ 
bles smaller than 0.01 mm^ can’t be discriminated from 
labeling artifacts m, and thus, are discarded. 



FIG. 2. 3D volume representation of two instances in the 
foam flow. The PBs and vertices are colored in yellow and 
blue for time steps to and to ± 35 s respectively. The scale bar 
is 300 pm. Red arrows indicate the flow direction. 


Velocity field - From the bubble tracking, we could 
measure their velocity around the obstacle. Statistics 
are performed in the diametral {rz) plane of the cylin¬ 
drical coordinates. We have checked that our results 
do not depend significantly on the angular coordinate 
0 (see also below), and we have averaged over this coor¬ 
dinate, as well as over time, thanks to the steadiness of 
the flow. The diametral plane is meshed into rectangular 
boxes (0.25 x 0.40 mm^). Gonsistently with the angular 
averaging procedure, we have checked that the number 
of bubbles analyzed per box is roughly proportional to 
the distance of the box to the symmetry axis (data not 
shown). Averages are weighted by the bubble volumes. 
The velocity field is plotted in Fig. In average over all 
patches, the 6>-component of the averaged velocity vector 
is 50 times smaller than its rz-component, hence the flow 
is axisymmetric. As expected, the velocity is uniform far 
from the obstacle, its amplitude decreases close to the 
leading and trailing points of the sphere, and increases 
along its sides. Accordingly, there is a clear recirculation 
zone surrounding the obstacle in the frame of the flowing 
foam (Fig. [^). It is worth noting that, compared to 2D 
foam flows around a circular obstacle [iBHin], the range 
of influence of the obstacle on the flow field is smaller. 

Interestingly, there is a zone downstream the obsta¬ 
cle and close to the symmetry axis where the streamwise 
velocity component is larger than the mean velocity or, 
equivalently, where the velocity opposes that of the ob¬ 
stacle in the frame of the flowing foam. This reminds 
the so-called negative wake^ revealed in viscoelastic flu¬ 
ids m and also evidenced in 2D foams m- However, 
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FIG. 3. Velocity fields in the (rz) plane, (a) in the lab 
frame. The normalized velocity held obtained by subtracting 
the mean how velocity is shown in (b). The half-circle repre¬ 
sents the obstacle (diameter 1.5 mm). The red arrow centered 
on the semi-obstacle gives the velocity scale of 8 /im/s. Blue 
arrows show the negative wake effect. 
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FIG. 4. Projection of the bubble deformation held in the (rz) 
plane. Ellipses of bubbles are dilated by a factor of 3. The 
colormap gives the amplitude of the normalized deformation 
in the azimuthal direction. 
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a difficulty intrinsic to the 3D axisymmetric geometry is 
that the statistics is poor in these boxes close to the sym¬ 
metry axis (about 10 bubbles per box over the full run), 
and should be improved in the future. 

Bubble deformation - Given the set of coordinates {r} 
of the voxels inside a bubble, we dehne its inertia ten¬ 
sor I = ((r — (r)) (8) (r — (r))), and its shape tensor as 
S = This operation is valid because I (and hence 

S) is a symmetric and dehnite tensor. The bubbles de¬ 
formation is quantihed by the eigenvectors/values of the 
shape tensor. In good approximation, two of them (S'/^ 
and are found inside the (rz) plane, the other cor¬ 
responds to the projection of the tensor along the az¬ 
imuthal direction (So). An effective radius is defined by 
^eff = (Sf'^Sf^So)^^^. The bubbles deformation in the 
(rz) plane is represented by ellipses of semi-axes 5'/^ and 
S~^. The direction of the largest one, 5'/^, is emphasized 
by a line across the ellipse (Fig. [^. Deformation in the 
azimuthal direction is quantified by (Sq — in 

colormap. The orientation of the ellipses in the (rz) plane 
exhibits a clear trend comparable to the 2D case [niiiT]. 
They are elongated streamwise on the obstacle side and 
at the trailing edge. In between, the ellipses rotate 180° 
to connect these two regions. We noticed that the defor¬ 
mation of the bubbles is much smaller than for a 2D foam 
with the same liquid fraction [min]- The deformation 
in the azimuthal direction exhibits dilation/compression 
up to 10% only. The quantity (Sq — is posi¬ 

tive upstream (oblate shape) to favor the passage around 
the obstacle (Fig. [^. The third dimension tends there¬ 


fore to reduce the bubble deformation in the (rz) plane 
by increasing the deformation in the azimuthal direction. 
This effect is opposite downstream, right after the obsta¬ 
cle, where the bubbles are prolate. 

Plastie rearrangements - Automated tracking of bub¬ 
ble rearrangements was hindered by the high sensitivity 
of such procedure to small defects in the reconstruction of 
the bubble topology. Description of the contact between 
bubbles requires to rebuild precisely the faces between 
bubbles, which would require a finer analysis [22]. Nev¬ 
ertheless, we managed to detect manually four individual 
events, corresponding to the rearrangements of neighbor¬ 
ing bubbles. We provide below a detailed description of 
one typical example (Fig. |^; the features of the three 
other ones were found to be the same. Those rearrange¬ 
ments consist of the swapping of four neighboring bub¬ 
bles, with an exchange of four-sided faces, called Tls or 
quadrilateral-quadrilateral (QQ) transitions by Reinelt 
and Kraynik [23l [24| . We did not observed three-sided 
faces during a T1 as reported by [25]. These are likely 
highly unstable, transient states which are too short-lived 
to be captured by tomography. The QQ transitions ob¬ 
served involve two bubbles losing one face and two bub¬ 
bles gaining one face. As can be seen on the projection 
plane across these four bubbles (Fig. and d), this is 
analogous to Tls in 2D, which always involve four bub¬ 
bles, two losing one side and two gaining one side. The 
distance between the two bubbles coming into contact 
decreases of 150 /im, from 1.10 mm before the T1 to 
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FIG. 5. Example of a T1 in 3D. The white and cyan bubbles 
lose contact, whereas the green and yellow bubbles come into 
contact. The red bubbles are the two bubbles that are in con¬ 
tact with these four bubbles over the rearrangement. Three 
different projections are shown: across the four bubbles that 
swap neighbors, (a) before and (d) after the Tl; (b) in the 
plane of the face that is about to disappear, and (e) in this 
plane after the Tl; (c) in the plane of the face that is about 
to appear and (f) in this plane after the Tl. 


0.95 mm after, while the distance between the two bub¬ 
bles losing contact increases of 200 /rm, from 1.03 mm 
before the Tl to 1.23 mm after. The distance between 
the two closest bubbles to the four ones directly involved 
in the Tl changes much less, from 1.48 to 1.52 mm. This 
corroborates the vision of a Tl quite similar as in 2D, 
acting as a quadrupole in displacement, with most effect 
on the bubbles in the plane. On the other hand, the vari¬ 
ation of shape anisotropy of the bubbles involved in the 
Tl did not show significant trends. 

We went further on in the characterization of the spa¬ 
tial structure of those rearrangements. Bubble faces com¬ 
prise a thin film surrounded by a thick network of PBs 
and vertices. All the faces that we observed are close 
to be isotropic, i.e. they do not show any significant 
elongation in a given direction. We have observed that 
the thin film part is usually very small for faces that 
are about to disappear, or that have just appeared, dur¬ 
ing a rearrangement. On the other hand, due to the 
finite radius of the PBs and of the finite size of the ver¬ 
tices, the “skeleton” of these faces is not arbitrarily small. 
Quantitatively, we measured on the images a PB radius 
i?c = 130 fim. We also measured the area of the skeleton 
of the faces on 2D projections along the plane of the faces 
(we did not observe significantly non-planar faces). We 
always found skeleton areas larger than 3.4 x 10^ /rm^, 
which is of the order (a bit larger) of R^. This suggests 
an interesting analogy with the cut-off edge length in 2D 
foams expected in theory [261127], and measured in both 
simulations [28] and experiments [16]. In 2D foams and 
emulsions, when the distance between two approaching 


vertices reaches a certain length, a rearrangement occurs. 
This happens usually when the two PBs decorating the 
two neighboring vertices start to merge; hence, the order 
of magnitude of the cut-off length is Rc. For 3D foams, 
our observations suggest that there is a cut-off area of 
the order of R^ below which a face becomes unstable, 
triggering a rearrangement. 

In summary, we have provided the first experimental 
measurement of a 3D time- and space-resolved foam flow 
measured directly from individual bubble tracking, with 
novel results on all the essential features of liquid foam 
mechanics: elasticity, plasticity and flow, through de¬ 
scriptions of shape field, Tl events, and velocity field. 
Such experimental results could be achieved thanks to 
the recent advances of both high resolution and fast X-ray 
tomography and quantitative analysis tools. Perspectives 
include further refinements of the reconstruction meth¬ 
ods OEa, to fully automatize the detection of rear¬ 
rangements, to increase statistics and to study various 
geometries. Imaging the 3D flow at the bubble scale may 
shed new light on pending issues on shear localization [9] 
and nonlocal rheology [29]. 
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